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ABSTRACT 

In this work we analyze the physical properties of a sample of 153 star forming galaxies at z^0.84, 
selected by their Ha flux with a NB filter. B-band luminosities of the objects are higher than those of 
local star forming galaxies. Most of the galaxies are located in the blue cloud, though some objects are 
detected in the green valley and in the red sequence. After the extinction correction is applied virtually 
all these red galaxies move to the blue sequence, unveiling their dusty nature. A check on the extinction 
law reveals that the typical extinction law for local starbursts is well suited for our sample but with 
E(B-V) s t ars =0.55 E(B-V) 9as . We compare star formation rates (SFR) measured with different tracers 
(Ha, FUV and IR) finding that they agree within a factor of three after extinction correction. We find 
a correlation between the ratios SFRfc/v"/SFR# q , SFR/fl/SFR^Q, and the EW(Ha) (i.e. weighted 
age) which accounts for part of the scatter. We obtain stellar mass estimations fitting templates to 
multi- wavelength photometry. The typical stellar mass of a galaxy within our sample is ^10 10 M Q . 
The SFR is correlated with stellar mass and the specific star formation rate (sSFR) decreases with it, 
indicating that massive galaxies are less affected by star formation processes than less massive ones. 
This result is consistent with the downsizing scenario. To quantify this downsizing we estimated the 
quenching mass Mq for our sample at z^0.84, finding that it declines from Mq ~10 12 M q at z~0.84 
to Mq ~8x10 10 M q at the local Universe. 
Subject headings: galaxies: evolution - galaxies: high redshift - galaxies: starburst 



1. INTRODUCTION 

During the last two decades the cosmic Star For- 
mation History of the Universe has been widely stud- 
ied in order to better constraint galaxy formation 
and evolution models. Large-area surveys as well 
as the use of larger telescopes have consolidated our 
kno wledge at low-inte rmediate redshifts (z=0. 0-1.0) 
(see [Hopkins & Beacom 2006). Several measurements 
exist at higher redshifts dPerez-Gonzalez et al.l 120051: 



exist at higher rcdshitts (Fcrcz-Lronzalcz ct al. 
Geach et all [20081: iReddv fc Steidell l2009t lHaves et al.l 
201C0 and at present we ha ye started to probe the mo st 



distant Universe at z~7-8 (jBouwens et al.ll2009l 12010). 
In general, the Star Formation Rate density (SFRd) 



history from the Local Universe to z^l is well accepted. 
Near z~l, where the rise in SFRd from the local Universe 
slows down, there are several measurements obtained 
through samples selected in a variety of ways and using 
different Star Formation Rate (SFR) tracers. Results 
from Ha, UV and IR agree reasonably, alt hough with 
highe r dispersion than at lower redshifts (|Garn et al.l 
2010) . This scattering is originated (at least partially) 
due to the two aforementioned processes: the sample se- 
lection or/and the SFR estimation. The estimations of 
the SFRd at t his redshift have been measured mainly 
through UV (ILillv et all fl996t IConnollv et al.l fl997l: 
Cowie et al.lll999HWilson et al.ll20 02: Schimin ovich et all 
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and Ha ( Glazcbroo k et al 
1999; 'Yan et alj|1999t IHopkins et al.l[2000t jTresse eTal 
2002; Do hertv et al.|2006HVillar et alil200cl [Sobral et a! 
2009; Lv et al.ll201lD . Therefore, it is important to con- 



straint the potential differences that arise when one or 
another tracer is used to estimate SFRs. 

Estimations of SFR through the UV flux are very sen- 
sitive to the extinction correction. The UV also probe 
older populations than Ha, thu s being more sensit ive to 
recent star formation history (|Calzetti et alJ l2005). The 
IR on the other hand is not affected by exti nction but 
is very model dependent (|Barro et al.l 120 11 al ibi). More- 
over, it is not well understood how the old population 
of stars contributes t o the IR emission, but it could be 
a sign ificant fraction (|da Cunha et al.l 1200 8: Sali m et al.l 
2009). The Ha line is one of the best estimators as it 
is sensitive only to very young stars, not affected by re- 
cent star formation history, which may still be detected 
by UV or IR. It has the problem that at z>0.5 it moves 
to the near Infrared (nIR) domain, where large amounts 
of spectroscopy are still difficult to obtain, though sev- 
eral nIR multi-object spectrographs for 8-10 meter class 
telescopes are coming in the next years. In addition, 
different methods have been used for measuring the to- 
tal Ha line flux, which could lead to discrepancies if 
some effects are not properly corrected. On one hand we 
have spectroscopy, long- or multi- slit or through fibers, 
where aperture correctio ns are needed to recover the to- 
tal fl ux (see for example iDohertv et al.|[2006t lErb et ahl 
l2006f) . On the other hand we have the slit less spec- 
troscopy (lYan et al.lll99 9; Hopk ins et al.l 
row band (Villar et al. 2008; Sobral et al. 



2000) and nar- 
20091 ILvet al.l 
1201 If ) techniques, which have the advantage that the to- 
tal flux of the object is recovered and no aperture cor- 
rections are needed. Although it is needed to correct 
from the Nitrogen contribution if the filter is not narrow 
enough, this effect is easier to estimate than aperture 
corrections. Obviously, the extinction is also important 
at this wavelength, though less than in the UV. 

Thus, the Ha estimator is the best suited to study the 
instantaneous star formation. As mentioned before the 
problem is that at this redshift the line is observed in 
the nIR and little data is available up to date. Then, it 
is also interesting to asses if UV and IR provide SFRs 
comparable to those of Ha, given the large amount of 
data in these wavelengths available today. 

The selection methods of the samples are also different 
and target different populations. Selecting the sample in 
the UV for example, implies a bias against very obscured 
galaxies, which may not be detected unless very deep 
observations are carried out. On the other hand, the IR 
selected samples will favor the selection of objects with 
large amounts of dust, missing the blue and dust-free 
objects. Selections based in the Ha line will select the 
objects with ongoing star formation and thus only star 
forming galaxies (and AGN) will be selected. 

Thus, to study the population of galaxies that are ac- 
tively forming stars it is necessary to have a well defined 
sample of star forming galaxies and one of the best tech- 
nique today is the use of narrow band filters targeting 
Ha. 

A population of star forming galaxies selected in this 
way is ideal to study the SFR sequence, which refers 
to the correlation that exists between SFR and stellar 
mass. This correlation has been found at a wide range of 
redshifts although there has b een found evolution with 
respect to this parameter (see lDutton et al.ll2010l for a 
review). While the slope is almost constant, the SFR ze- 



ropoint increase from the local Universe to redshift z^2. 
From this redshift on, the trend remains almost constant 
with little evolution. There exist some discrepancy in the 
slope, which is s omewhat lower than unity in all cases 
(see for example INoeske et al.l 1200711 lElbaz et all 120071; 
ISalimet ai.ll2007ft . 

However, some works do not find this correlation. 
iCaputiet all (|2006f ) do not find any trend fo r their MIPS 
selecte d sample at z~2. More recently, iSobral et al.l 
(2010) did not find any evidence of this sequence for their 
HiZELS sample at z~0.84, selected with a narrow-band 
filter targeting Ha. 

This correlation implies that the slope of specific SFR 
(sSFR) versus mass is higher than -1. Obviously, if no 
correlation between stellar mass and SFR is found the 
slope of sSFR vs. stellar mass is simply -1. The slope of 
this correlation is very important, as it tells us how de- 
creases the importance of star formation over the already 
formed stellar mass. 

In this work we use a narrow band select ed sample of 
star fo rming galaxies at z~0.84, presented in lVillar et al.l 
(2008), to compare SFRs obtained from different tracers 
and to study the relation between stellar mass and SFR. 
The sample is very well suited for this study as: i) it is 
directly selected by star formation, so we are not biasing 
the population of star forming galaxies; ii) the use of 
the narrow band filter technique provides us reliable Ha 
SFRs to compare with FUV and IR estimations. 

This paper is structured as follows. In Section [21 we 
present the sample and the available datasets. In Sec- 
tion [31 we describe the methods used to get rid of AGN 
contaminants. Absolute magnitudes and color are pre- 
sented in Section |4j Section [5] presents a comparison of 
SFRs obtained through different estimators as well as a 
check on the extinction law more suited to our sample. 
Section [5] presents the stellar masses and their relation 
with star formation rate. Finally, we summarize our re- 
sults and conclusions in Section [71 

Throughout this paper we use AB magnitudes. We 
adopt the cosmology Hq — 70 km s _1 Mpc 1 , ft m = 0.3 
and fl\ = 0.7. 

2. DATA 

2.1. Sample 

This paper analyzes an Ha selected sample of galaxies 
at z=0.84. The objects are selected by their emission in 
the Ha+[NII] line and are thus selected due to intense 
star formation (except when activity at n uclei level is 
resen t). The sample was first described in lVillar et al.l 
20081 hereafter V08) and the reader is referred to that 
paper for full details on the sample selection criteria. A 
brief summary of the process is presented here. 

The sample was built using narrow and broad band 
images in the J band of the near infrared. The narrow- 
band filter used in this work is J-continuum ( Jc) centered 
at 1.20/im, corresponding to Ha at z=0.84. The search 
was performed using the near-infrared camera OMEGA- 
2000J of the 3.5m telescope at the Calar Alto Spanish- 
German Astronomical Center (CAHA). OMEGA-2000 is 
equipped with a 2k x 2k Hawaii-2 detector with 18/im pix- 
els (0"45 on the sky, 15'xl5' field of view). Three point- 
ings were observed, two in the Extended Groth Strip 

1 http://www.mpia-hd.mpg.de/IRCAM/O2000/index.html 



(EGS) and another one in the GOODS-North field, cov- 
ering a whole area of ^0.174 deg 2 , and reaching 70% 



1 cm" 2 



completeness at a line flux of ~1.5xl0 erg s 

In a first step, 239 emission line candidates were se- 
lected (once excluded the stars) by their flux excess in 
the narrow band, showing a J-Jc colour excess signif- 
icance n CT > 2.5 in one or several apertures. Spectro- 
scopic and photometric redshifts were then used to rule 
out contaminants, either emission line galaxies at other 
redshifts or objects selected by spectral features or noise. 
First, the sample was cross-checked against redshift cata- 
logs on GOODS-N and EGS fields. A total of 76 objects 
were confirmed as genuine Ha emitters in the narrow 
band redshift range, 43 in the Extended Groth Strip and 
33 in GOODS-N. Contaminants were mainly emission 
line galaxies at other redshifts, including a small sample 
of [OIII]AA4959,5007 emitters at z~1.4, and a few ob- 
jects not selected by line emission. The accuracy select- 
ing emission line galaxies was very high, around ~90%. 
Spectroscopic redshifts were only available for 98 objects 
in the sample, therefore photometric redshifts were used 
to get rid of interlopers for the rest of the sample. The 
quality of the estimated photometric redshifts was very 
high, with 86% of the objects in the whole EGS and 90% 
of the objects in GOODS-N (with reliable spectroscopic 
redshifts) within cr z /(l + z) <0.1. Considering the pho- 
tometric redshifts, a total of 89 objects, 64 in the EGS 
and 25 in GOODS-N, were added to the final sample. 

Since the original paper was published new spectro- 
scopic data has been made public, increasing the num- 
ber of confirmed sources in 18 objects, for a total of 94. 
Only two objects are found to be incorrectly classified as 
Hq emitters, which has been removed from the sample. 
Thus, the selection efficiency found in the original sample 
remains very similar. 

The final sample of Ha emitters at z=0.84 contains 165 
objects, 107 in the EGS and 58 in GOODS-N, 94 (57 
%) of them confirmed by optical spectroscopy (after in- 
cluding three objects with low quality spectroscopic red- 
shift). However, due to insufficient complementary data, 
we have discarded 12 objects. Hence, the sample used in 
this paper is composed of 153 objects. 

Line fluxes ha ve been recom p uted using the formal- 
ism described in IPascual et al.l (J2007D (hereafter P07), 
introducing the Nitrogen contribution in the filters ef- 
fective widths. This forces us to assume an initial value 
for the nitroge n con tribution, setting it to the average 
value found in |pl: J([iVi7]A6584)//(#a)=0.26. This 
provides an initial estimate of the Ha line flux, without 
nitrogen contribution. With the equivalent width, we can 
estimate the nitrogen contribution, given the correlation 
between EW(Ha) and I {[N I I]X658A) / 1 (H a) found in 
the local SDSS sample (see P07). We then reestimate 
the nitrogen contribution and compute again the line 
flux and equivalent width. The latter provides a new 
estimation of the Nitrogen contribution. The process is 
repeated until it converges, usually in two or three steps. 

The line flux estimation through the narrow and broad 
band filters assumes a most likely redshift for the object, 
i.e. we assume a redshifted wavelength for the emission 
line based in the shape o f the filter and in the cosmology 
(see IPascual et al.1 l2007V In fact, the objects distribute 
along the wavelength range covered by the narrow band 
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Fig. 1. — Line flux ratio between those computed with the indi- 
vidual spectroscopic redshift and with the average redshift. The 
narrow band filter shape is represented by the continuous line. The 
shaded region comprises the redshift range where the filter trans- 
mission is above 80%, where most of the objects are detected. The 
dotted line represents the Ha line surrounded by the Nitrogen lines. 



filter, being most likely detected near the filter's cen- 
tral wavelength, where the transmission is high. Thus, 
assuming a most likely redshift for the objects seems rea- 
sonably for the majority of them. However, for the ob- 
jects that are selected in the wings of the filter, where the 
transmission falls abruptly, the recovered fluxes differ sig- 
nificantly. Even in the regions of high transmission, there 
could be important effects due to the presence of the ni- 
trogen lines. We can correct these effects introducing 
the real redshift, which we know for half the sample, in 
the equations to compute the line flux (see IPascual et all 
12001 . 

In figure [T] we compare the line fluxes estimated with 
the real redshift versus the ones estimated in the gen- 
eral way. The narrow-band filter shape is also shown as 
reference. It can be seen that, for the objects that fall 
in the wings of the filter, the line flux is clearly subes- 
timated. It is worth noting also that a significant line 
flux fraction is lost even when the transmission is high, 
as is the case of the objects with higher redshifts within 
the shaded region in the figure, where the transmission 
is always above 80%. This is due to the fact that the 
[NII]A6584A line, which is the most intense of the two 
nitrogen lines, shifts to wavelengths where the transmis- 
sion of the filter is very reduced, and hence, if we do not 
take this effect into account, we overcorrect the nitrogen 
contribution, estimating fainter line fluxes. The effect is 
also present at shorter wavelengths, though in that case is 
the other Nitrogen line ([NII]A6548A) the one shifted to 
wavelengths with lower transmission. However, as this 
line is 3x weaker than the other one, the effect is less 
pronounced. 

The amount of extinction for each galaxy was esti- 
mated through the FIR to UV flux ratio or the UV slope 
when the FIR data was not ava ilable (see V08). W e 
used the extinction law derived bv lCalzetti et al.l (|2000f) . 
As new data is available now, specially regarding MIPS 
24/xm, we have recomputed the extinctions for all the ob- 
jects, considering, in addition, the results on the check 



on the extinction law (see section f5 . 1 1) . The median ex- 
tinction for our sample is 1.24 m in Ha, adopting values 
between m and 3.8 m . Once applied these corrections, 
Ha luminosities and star formation rates were computed. 

2.2. Additional data 

In order to estimate the different properties ana- 
lyzed in this paper, we use several additional data sets 
sampling a wide range of the electromagnetic spec- 
trum, from the Far Ultraviolet (GALEX FUV) to the 
Mid Infrared (MIPS 24/im). These complementary 
data sets have been colle cted as part of the UCM 
Rainbow database (see V08. lPerez-Gonzalez et al1l2008t 
IBarro et al.H2011al for details), and have been gathered 
in part by the AEGIS (IDavis et alJl200l and GOODS 
(jDickinson fc GOODS Legacy Team! 120011 ) projects. 

Briefly, in the EGS we have used optical data gath- 
ered with MegaCam at the 4m Canada France Hawaii 
Telescope (CFHT), covering the following bands: u* , g', 
r', i\ and z'. We have also used B, R, and I images 
obtained with the same tel escope but with the camera 
CFHT 12k, as described in lCoil et al.l (|200l . Deep op- 
tical R band data taken with SuprimeCam at Subaru 
8m as part of the Subaru Su prime-Cam Weak-Lensing 
Survey ([Mivazaki et al.l [20071 ) are also available. In the 
domain of the near infrared, images in the J band were 
obtained with Omega2000, as part of the data necessary 
to make the sample select ion, and K-band images ob- 
tained with Omega prime ([Barro et al.l [2009). Both in- 
struments were located at the 3.5m telescope at CAHA. 
Space-based optical images acquired with the Advanced 
Camera for Surveys (ACS) on board HSTaze available in 
two bands: F606 W and F814W (hereafter Vso 6 and isu). 
The Galaxy Evolution Explorer (GALEX, iMartin et al.l 
|2005[ ) provides ultraviolet deep images in the far ultravi- 
olet (FUV; 153 nm) and the near ultraviolet (NUV; 231 
nm). The space observatory Spitzer observed the EGS 
field at 3.6, 4.5, 5.8, and 8 j Ltm with the IRAC in strument 
and in 24 //m with MIPS (jBarmbv et al.l 120081 ). 

In GOODS-N we made use of deep optical and near 
infrared images (UBVRIzHK s , iCapak et al1l2004l ). as 
well as our own J- and K-band images , both of them ob- 
tained with Omega2000 (see V08 and JBarroet al.l 120091 
for details). As in the case of the EGS, space based obser- 
vatories provide us with ultraviolet and infrared data, as 
well as high resolution additional optical data. GALEX 
observed the region in Far- and Near- Ultraviolet chan- 
nels. Spitzer observed in the mid infrared (3.6 to 8 ^m; 
IRAC) and in the far infrared (24 ^m; MIPS). The ACS 
on board the HST adds optical data in four bands F435B 
(£435), F606W(V 606 ), F775W(i 775 ), and F850LP (z 850 ). 

3. AGN CONTAMINANTS 

The selection of a sample through the Ha line is sensi- 
ble to Active Galactic Nuclei (AGN) contamination, as 
they are also powerful emitters in this line. Although 
both AGN and SFR could contribute together to the 
flux, disentangling both components is out of the scope 
with the available data. Thus, we will remove the objects 
classified as AGN from our sample. 

In this work we detect the presence of 13 (8%) AGN us- 
ing two complementary methods: X-ray luminosity and 
mid-IR colors. 
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Fig. 2. — IRAC color-color plots for the selected sample with 
photometry in the four bands. The wedge delimited by the dashed 
polygon encloses the emitters powered by an AGN (open triangles). 
Seven objects (filled triangles), that fall outside this wedge, with 
positive [3.6]-[4.5] color have been also considered AGNs because 
they have a rising SED and the errors make them compatible with 
being located inside the wedge. 

3.1. X-ray luminosity 

We have cross-correlated our sample with the avail- 
able X-ray catalogs in the EGS and GOODS-N fields. 
In the EGS fields we have the AEGIS-X X-ray catalog 
([Laird et al.ll2009l ). which covers a large area within the 
EGS. Observations were made with the Chandra X-ray 
observatory with nominal exposure time of 200ks. In 
GOODS-N we have used the catalog created by Laird et 
al. from observat ions taken by Chandra with an expo- 
sure time of 2Ms ([Alexander et al.l 120031 ). We find three 
X-ray counterparts in the EGS and four in GOODS-N, 
within a 2" search radius. The three objects in the EGS 
present high X-ray fluxes (Lx >6xl0 42 erg s _1 ), reveal- 
ing their AGN nature. In GOODS-N, due to the depth 
of the observations, we find three objects whose X-ray 
luminosities are compatible with an star formation ori- 
gin. The der i ved S FRs, using the calibration given by 
iRanalli et al.1 (J2003), agree within a factor of three with 
the Ha derived ones. Thus we have only discarded the 
four objects with X-ray luminosities whose origin could 
only be attributed to an AGN. 

3.2. Mid Infrared colors 

Some AGN are heavily obscured and even the deep- 
est X- ray observations can miss a significant fraction of 
them (Park ct al. 2010). In this case the X-ray emission 
is absorbed by the circumnuclear dust and re-emitted in 
the infrared. One way to detect these obscured AGNs 
is lookin g at their m i d-IR colors. The color criterion de- 
fined bv lStern et al.l ([2005D is based in the differences in 
the mid IR emission shown by star-forming galaxies and 
AGNs. The star-forming gala xies SED peaks a t 1.6^m, 
falling at longer wavelengths ([Garn et al.l 120101 ). In the 
case of an AGN the emission do not decreases at longer 
wavelengths, due to the re-emission of light absorbed by 
the circumnuclear region in the mid-IR. Unfortunately 
the distinction be comes less p r onoun ced at our redshift, 
as pointed out bv lStern et al.l (|2005|) . 

In figure [2] we show the objects that fulfill the crite- 
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Fig. 3. — Histogram of rest-frame B-band absolute magnitudes. 
The thick line shows the distribution for the z~0.84 sample. The 
dashed line correspond to the UCM local sample, selected also by 
their Ha emission. 

rion (inside the dashed polygon), which are AGNs, as 
well as the rest of the objects, which are pure star form- 
ing galaxies. A t otal of thre e galaxies fall within the 
wedge defined by I Stern et all Seven other objects with 
positive [3.6]-[4.5] color fall relatively close to this wedge 
(except one), although do not fulfill the criterion. We 
have decided to consider these objects as AGN contam- 
inants, given that photometry errors could have placed 
them outside the region and that they have a rising SED 
([3.6]-[4.5]>0). Only one of these ten objects have an 
X-ray counterpart. 

We note that this classification could be selecting star 
forming galax ies in stead of AGNs, as pointed out by 
iDonlev et al.l (|2008l ). We decided to exclude them all 
to make sure we are not introducing any AGN. 

Another way to check the prese nce of obscured AGN 
is thr ough the power-law criterion (|Alonso-Herrero et al.l 
I2006T) . We do not find any galaxy showing this charac- 
teristic power-law shape in the mid-IR SED. 

4. PHOTOMETRIC PROPERTIES 

The histogram of rest-frame B-band magnitudes for 
our sample is shown in figure [3) The median of the 
distribution is Ms=-20.5 m , reaching the most luminous 
objects M#=-22.5 m . The standard deviation of the dis- 
tribution is 0.9 m . For comparison we also show in the 
figure the UCM local sampl e of star forming galaxies 
iZamorano et al.l (|1994 I1996T) . selected also by the Ha 
line flux. In the local sample no galaxies brighter than 
Mb ~-22 are detected, suggesting that star forming 
galaxies at z^0.84 are more luminous than their local 
analogous. The volume sample in each survey is very 
different: 10 5 Mpc 3 for the typical object in the UCM, 
while for our sample the surveyed volume is ~15xl0 3 
Mpc 3 . However, given the larger volume explored in the 
UCM survey with respect to our survey, it is clear that 
star forming galaxies at z^0.84 are in general brighter in 
the B band. 

A clea r bimodality in the c olor of galaxies was first 
found bv lStrateva et al.l (|200lD analyzing the optical col- 
ors of the SDSS sample. Galaxies divide mainly in 
two groups: the blue cloud and the red sequence. The 
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Fig. 4. — Histogram of rest-frame NUV-R colors for our sample. 
The black line represents the sample without applying the correc- 
tion for extinction. The dashed line shows the distribution once 
the extinction has been corrected. It can be seen that most ob- 
jects fall in the blue cloud ((NUV-R)<3.5). Before the extinction 
correction is applied some objects fall in the green valley (3.5< 
(NUV-R)<4.5) and in the red sequence ((NUV-R)>4.5). After it 
is applied only two objects remain outside the blue cloud. 

blue cloud is populated by star forming galaxies whereas 
galaxies with no recent star formation fill the red se- 
quenc e. This b imodali ty i s also present at higher red- 
shifts (jWillmer et al.ll2006t iFaber et al.ll2007l ). An inter- 
mediate region , the green valley, w as identified using the 
NUV-R color (|Wvder et all 120071) . The galaxies within 
this group are either in a transition phase from the blue 
cloud to the red sequence, due to the shutdown of star 
forma tion: or are star forming galaxies with high extinc- 
tion (jMartin et alll2007t ISaTim et all 120071 ). In figure H 
we depict the NUV-R color for our sample. Most of the 
sample belong to the blue cloud (NUV-R<3.5), in agree- 
ment with their star-forming nature. However, some 
galaxies fall in the green-valley (3.5<NUV-R<4.5) and 
a few of them in the red sequence (NUV-R>4.5). This 
can be explained due to the high extinction present in 
these galaxies and, indeed, when the extinction is cor- 
rected only two galaxies fall outside the blue cloud. One 
of them is not confirmed by optical spectroscopy so it 
might not be a real z~0.84 emitter. In fact, its photo- 
z x 2 distribution does not present a clear peak at that 
redshift but a flatter distribution. The other one, al- 
though confirmed by optical spectroscopy, is very close 
to another galaxy (<2") and its photometry might be 
affected. 

5. STAR FORMATION RATES 

This work uses the Ha luminosity as the principal es- 
timator of the instantaneous SFR of galaxies. The Ha 
line flux has been used to select the sample, making it 
very suitable to study star formation processes, as it has 
been selected by this property. However, star formation 
involves physical processes whose imprint become ob- 
servable along a wide range of the electromagnetic spec- 
trum: X-rays, ultraviolet, forbidden recombination lines 
([O II]A3727), far infrared, radio, etc. In this work, in 
addition to Ha, given the depth and coverage of the multi 
wavelength data available, we estimate SFRs through Far 



G 



Ultraviolet and Far Infrared luminosities. Each tracer is 
affected by different phenomenons and is originated by 
different physical mechanisms, related (at least in part) 
with star formation processes. Thus, the different results 
obtained with different tracers could yield some informa- 
tion about the properties of the galaxy that hosts the 
star formation processes. 

The Ha line is produced due to the recombination pro- 
cesses in ionized hydrogen present in the clouds of gas 
and dust that surround the newly formed stars. The 
massive stars of type O and B produce an intense radi- 
ation field capable of ionize the hydrogen atoms. When 
the equilibrium is reached, the recombination of the free 
electrons with the ionized hydrogen produce several emis- 
sion lines, being the Ha line one of the most luminous 
in the visible. To obtain the star formation rate from 
the Ha luminosit y (L# a ) we apply the relation given by 
lKennicuttl(|1998t l: 

SFR Ha (M Q yr- 1 ) = 7.9 x l(n 42 L ff a (erg s~ l ) (1) 

where a iSalpeteil (1955) IMF has been considered. Ha 
traces directly the SFR and has very low dependence on 
metallicity or on ionization conditions of the gas cloud. 
Among the adverse effects the most important are the 
extinction and the escape fraction of ionizing photons. 
The former is common to optical indicators and extinc- 
tions as high as ~4 magnitudes in the Ha line can be 
found in our sample, although the median extinction is 
1.24 magnitudes. The latter implies a subestimation of 
the star formation rate if the escape fraction is high. 
Fractions u p to 50% have been m easured for individual 
HII regions (|Oev fc Kennicut t 1997). However, this frac- 
tion turns out to be much lower when the whole galaxy 
is considered (whic h is our case) decreasin g to less than 
3% as measured by iLeitherer et al.l (119951) . Anot her ad- 
verse effect recently shown by iLee et all (|2009l ) is the 
underestimation of star formation rate for dwarf galax- 
ies. Nevertheless, this effect appears for star formation 
rates below 0.03 Mgjyr -1 , two orders of magnitude lower 
than our lowest star formation rate, so it does not affect 
the estimations for our samples. 

Ultraviolet emission comes directly from young massive 
stars formed in the star formation region. To compute 
the UV SFRs we use the following calibration iKennicuttl 
(fl998h : 

SFR FUV {M & yr- 1 ) = 1.4 x lQ- 28 L FUV (erg s^Hz' 1 ) 

(2) 
where Lpjjv is the FUV luminosity spectral density. Al- 
though we apply it to the UV flux in 1500 A, the calibra- 
tion is valid in the 1500-2800 A range, as the spectrum 
is nearly flat in that regime. 

The dust in a galaxy absorbs part of the radiation 
emitted at short wavelengths and re-emits it in the IR. 
This absorption is more intense at shorter wavelengths. 
Given that young stars radiate most of their luminos- 
ity in the ultraviolet, there exists a correlation between 
IR luminosity and star formation. The correlation be- 
tween luminous regions in Ha and in IR confirmed the 
validi ty of the latter as a valid star formation tracer 
(see iDevereux et al.|[T997l for details). More recently, 
observations carried out with the Spitzer space tele- 



scope, with improved resolution, confirmed those results 
(jCalzetti et alj|2005t iPerez-Gonzalez et al.ll2006ft . 

The calibration of the infrared emission as a star for- 
mation tracer is not simple as it depends on several fac- 
tors: geometrical distribution and optical thickness of 
the dust, fraction of emission coming from old stars, etc. 
In the ideal case the star forming regions would be sur- 
rounded by dust clouds, opaque enough to re-radiate all 
the region luminosity. However, these ideal conditions 
differ from the real scenario, due in part to the afore- 
mentioned factors. This complicates the calibration and 
increases the dispersion. 

Given the infrared luminosity L/fl(8-1000/im), the fol- 
lowi ng relation can be used to estimate the star formation 
rate (jKennicuttlll998l ): 

SFRjniMQyr- 1 ) = 1.71 x 10 10 L IR (L Q ) (3) 

where L/^ corresponds to the total infrared luminosity 
between 8 and 1000 /im. 

The great advantage of this tracer is that it is not af- 
fected by extinction. However, the estimations are sen- 
sitive to other factors like the dust spatial distribution, 
old stellar population contribution, etc. Moreover, it is 
necessary to estimate the total infrared luminosity in 
the range 8-1000 /im. In general, this is done with a 
few measurements in the infrared, usually with wave- 
lengths shorter than 24/im. Recent works have demon- 
strated that a better relation exists between Ha and flux 
at 24/mi (see, for ex ample. IPere z-Gonz alez et al.l 2006: 
Calzetti et all 120071 : iRieke et al.l 120091 : iKennicutt et al.1 
2009f ). However, this wavelength is not available for our 
z=0.84 sample, because the MIPS observed 24/im data 
turns into restframe ~13/im. 

5.1. Dust attenuation 

In order to properly compare the SFRs it is necessary 
to correct the effect introduced by the extinction. In 
V08 we computed the extinction using the dust flux to 
UV flux ration and the UV slope when the IR data was 
not available. This provided us the attenuation in the 
FUV band, from which th e attenuation a t the Ha line 
was estimated assuming a iCalzetti et al.l (|2000l ) extinc- 
tion law, given the star-forming nature of our sample. 
This law considers that the nebular emission is more ex- 
tincted than the stellar emission: E(B-V) s t ars = 7 E(B- 
V) 9 as, with 7=0.44. However, this factor may be differ- 
ent for different populatio ns and /or d u st geo metry and 
may depend on redshift. iGarn et all (|2010l ) found for 
the S09 samp le that 7 ^0.5, slightly higher than the 
typical value. lErb et all (|2006l ) found that, in order to 
reconcile their SFR estimations in the UV and Ha for 
their z^2 UV selected sample, the same color excess 
has to be affecting bo th the gas and stars, i.e. 7 ~1. 
Yoshikawa et al. (2010), using a sample of BzK selected 
galaxies at z^2, found that their data was consistent 
with the original 7=0.44 value, although galaxies with 
low SFRs are consistent with 7=1. 

The 7 factor arises from the fact that the UV and the 
nebular emission have different spatial origins, due to the 
different population of stars each one is tracing. Whereas 
the nebular emission is originated from very massive and 
young stars, the UV emission is originated by less mas- 
sive and older stars. Thus, this factor might be different 



depending on the star formation history of the galaxies. 
M oreover, the ex t inctio n law may be different than that 
of iCalzetti et al.l ([2000). Thus, it is important to esti- 
mate this value for ou r sample and, if possible, to verify 
the suitability of the ICalzetti et al.l (with the same or 
different 7) extinction law for our sample. 

We tackle this problem estimating the extinction law 
in the UV regime and in the Ha line. Thanks to the 
large amount of optical broad-band data available we 
can obtain several estimations of the SFR (affected by 
extinction) at different wavelengths within the UV, in 
addition to the Ha estimation. If we assume that every 
different SFR estimation, once corrected for extinction, 
shall give the same SFR, we can obtain the extinction in 
each wavelength comparing to the total SFR: 



C p jjtotal np jDuncor -i n0.4 K{UV n ) E(B — V) a tar 



SFR 



total 



QT^TDuncor-ir\0.4:K,(Ha)E(B — V)ga 



(4) 
(5) 



where as SFR* otai we use the SFR given by the IR(8- 
1000/im), UV n represents each different UV wavelength. 
Thus, we can obtain k(A) for different wavelengths: 
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Calzetti et al. 2000 R-4.0 
Cardelliet al. 1989 R=3.1 
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Fig. 5. — Derived reddening curve for our sample (filled cir- 
cles). The continuos line represents the Calzetti extinction law 
with 7=0.55. The dotted line corresponds to the Cardelli extinc- 
tion law with 7=0.46. The open circles represent the reddening 
that would have the stellar continuum for each extinction law and 
their corresponding 7 values. 
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where we have everything related to the color excess in 
the gas E(B-V) gas through the 7 factor. At this point 
we are interested in the shape of the extinction law and 
is therefore necessary to get rid of the amount of extinc- 
tion, parametrized by E(B-V) ga s for each galaxy. We 
normalize then by the value at 6563A: 



/ SFR total 

1 & \ c p nuueor 

n 6563 (UV n )-- ^ Rtotal 



log I 



C Z7 ijuncor 



Kf>563{Ha) = l 



(9) 



(10) 



where 7 is not known and therefore we can only mea- 
sure empirically 7 • k(UV„). If we fit these values to a 
extinction law we can obtain this factor, which was our 
original goal. 

In our sample there are data available for four bands in 
the UV rest-frame within the range 1900-3000A, where 
the spectrum is almost flat on ce the dust effect has been 
removed and we can use the iKennicuttl (|1998[ ) calibra- 
tion. In both EGS and GOODS-N there are observa- 
tions at ^1950A and ^2400A. There exist an aditional 
third band in each field but with different wavelength: 
~2650A in EGS and -2950A in GOODS-N. Thus, we 
sample four different wavelengths in the UV. There are 
72 objects for which all the UV and infrared needed data 
is available. 



There are two things we want to check: a) which is 
the 7 factor a ppropriate for our sample? and b) is the 
ICalzetti et al.l extinction law suitable for our sample? To 
answer these questions we fit two different extinction laws 
to the data: the aforementio ned Calz etti et al.l extinction 
law and the iCardelli et all (1989) one. The fitting pro- 
cess gives us the 7 factor needed to make the UV and Ha 
measurements consistent and it allow also to check which 
extinction law provide a better agreement with the data, 
through the computation of the x 2 value of each fit. 

In figure[5]we show the 7 • «(UV„) obtained computing 
the median of the different UV measurements. We also 
show the Calzetti and Cardelli extinction laws that best 
fit to the data points. The Calzetti extinction law is more 
consistent with our measurements. We obtain the follow- 
ing x 2 values for each fit: Xc q;0 o= - 2 versus Xcar89= - 6 - 
These values are below one, due to t he large errors w e are 
working with. We consider that the ICalzetti et al.1 law is 
best suited for our sample, as the residuals are lower and 
we are dealing with the same uncertainties. In both cases 
a heavier attenuation in the nebular gas than in the stel- 
lar continuum is needed, i.e a 7 factor lower than one. 
We obtain 7 Cq ;oo=0.55±0.20 and 7 C ar89=0.46±0.17. If 
we do the analysis on a galaxy by galaxy basis we ob- 
tain similar results. The extinction analysis on individ- 
ual galaxies and the relation with other properties will 
be presented in a future paper. 

If we repeat this process discarding all objects not con- 
firmed by spectroscopy we obtain very similar results. 
In this case the number of galaxies is reduced to 57 
and we obtain 7c a ;oo=0-56±0.20 and 7car89=0.47±0.17. 
The x 2 values for each fit are now: Xc a ioo~^-^ versus 
X 2 7 arg9 =1.4. A comparison among results (from this and 
other sections) for the whole sample and that only con- 
taining spectroscopically confirmed objects is shown in 
Tabl<£D 

To summarize, the ICalzetti et aLl extinction law is well 
suited for our sample with 7=0.55, a value slightly higher 
than the original 0.44 value. We have assumed this ex- 
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Fig. 6. — Comparison of SFRs inferred from Ha luminosity and 
FUV luminosity. Top: No extinction correction applied. Bottom: 
Extinction correction applied to both tracers. Filled circles are 
objects confirmed by optical spectroscopy, whereas open circles 
are objects without spectroscopic confirmation. We show also the 
SFRFfyvVSFRj/ a ratio versus SFRj/a- 

tinction law with this 7 factor on the dust attenuation 
estimations. 

5.2. Comparison 0/ SFR tracers 
5.2.1. Ultraviolet 

In this section we compare the FUV derive SFR with 
that coming from the Ha luminosity. The reader should 
note that in this section we use the FUV luminosity 
(1500A), which has not been used in the computation 
of 7 in the previous subsection, thus assuring the inde- 
pendence of the results. 

The comparison between Ha and FUV star formation 
rates is plotted in figure [51 Top figure shows SFRs esti- 
mated without extinction corrections. Objects confirmed 
by optical spectroscopy are shown as filled circles whereas 



™^ 


_r 


_r 


: n. 



-2-10 1 2 

MSFRVSFRtfJ 

Fig. 7. — Histogram of the ratio SFRf(jv/SFR£f a in logarithmic 
scale. Most of the objects concentrate around the unity ratio. 
There are also wings on both sides of the distribution, with some 
extreme cases where SFR is over(sub)-estimated up to 10-100. The 
tail of objects for which the UV sub-estimates the SFR is more 
extended than in the case of the Ho. The dashed lines encloses the 
objects whose SFRs agree within a factor of three. 

the objects lacking spectroscopic confirmation are shown 
as empty circles. The effect of the reddening is clearly 
visible as FUV SFRs are, in general, lower than those 
obtained through Ha. The median value and standard 
deviation for the FUV estimates are (SFRjrj/v r )=l-5^o!9 
Moyr" 1 , while for the Ha line we find (SFR Ha )=3-5tf% 
MQyr -1 . Objects not confirmed by spectroscopy show a 
higher dispersion, although compatible with that of the 
confirmed objects once a few outliers are removed. 

In the bottom panel of figure |6j we show the effect 
of applying the extinction corrections. The SFR range 
spans considerably, going from 2-10 MQyr -1 when the 
effect of extinction is not corrected to 2-300 M Q yr _1 
when it is corrected. Estimations coming from both 
tracers now agree within a factor of 3. The statistical 
values are in this case: (SFR^uv^lOiy 1 Moyr -1 and 

(SFRHa)=ll-7 2 M Q yr~ 1 . The good agreement corrob- 
orates that our extinction corrections are working well 
and that these galaxies do not host star forming regions 
totally attenuated in the UV but visible in Ha, at least 
globally. There still can be regions totally obscured both 
in UV and in the optical, which will only arise in IR 
observations. We will explore this possibility in section 
ET2~2l 

To explore in more detail the differences between both 
tracers we study the SFR^ r/y/SFR# Q ratio for each 
object. The median value is (SFRF[/y/SFR# Q }=0.89, 
which tells us that the Ha line yields slightly higher val- 
ues than the FUV for the star formation rate, although 
compatible considering the errors. As we use the FUV 
luminosity, which has not been used in the computation 
of 7, it is possible to have ratios below or above one. As 
the FUV is at shorter wavelength, the higher extinction 
could totally attenuate more regions than at higher wave- 
lengths, thus underestimating the SFR. If we use SFRs 
obtained from 2800A instead, the ratio SFR28oo/SFR# Q 
becomes one, as this wavelength is in the regime used 
in the extinction law check. The number of objects is 



also different as objects used in the extinction section 
must have IR data. The distribution of ratios is shown 
in figure [7] Although the agreement is quite good, with 
90% of the objects within a factor of three, there are 
objects whose SFR is overestimated by Ha, with a few 
in the opposite case. If we consider only our spectro- 
scopic confirmed sample we find almost the same results: 
{SFR Ha )=Ut 2 9 3 Moyr^ 1 and (SFR FC/y /SFR HQ )-0.87. 
The general agreement b etween both trace rs is also 
found in the loc al Universe (jSalim et al.1 120071 ). At z~2 
lErb et all (2006) also compare d these tracers for a sam- 
ple of Lyman-break galaxies (|Steidel et al.l [T9 96. 1999) 
selected through the U n G R criterion (jAdelberger et al.l 
120041: iSteidel et all 12004(1 . Their result shows good 
agreement between both tracers with a dispersion sim- 
ilar to that of our samp le. At that same red- 
shift, lYoshikawa et al.l (J2010) find that both tracers are 
roughly consistent, although SFRs measured with Ha are 
systematically larger by 0.3 dex. 

5.2.2. Infrared 

Infrared SFRs are very interesting to check if Ha is 
losing substantial star formation due to dust attenuation. 
Top panel in figure [8] shows the comparison between both 
tracers before applying any extinction correction. Not 
surprisingly, Ha sub-estimates systematically the SFR. 
There is a large scattering, which reflects the different 
attenuation that suffer each galaxy. It is worth noting 
that we only have infrared luminosities for the fraction 
of the sample detected with MIPS (100/140, 71%). The 
IR limiting flux for a completeness of 80% is 83/iJy in 
our surveyed fields, which translates into ~10 M Q yr _1 
for our redshift. 

The median and standard deviation is 
(SFR H „}=4.1J:?^ Moyr- 1 for Ha while for the IR 

we obtain (SFRj fi )=15Jlfi M©yr _1 . As in the case of 
UV, we do not observe systematic differences between 
spectroscopically confirmed and not confirmed objects. 

Once the extinction corrections are applied (see bot- 
tom panel of figured]) we find (SFR# Q )=15±^ MqVt- 1 , 
which agrees very well with the IR derive value. If we 
work with the ratios of SFRs we find that Ha provides 
slightly higher estimates ((SFR/fl/SFR ffQ )=0.95), but 
in agreement with IR estimates within uncertainties. In 
figure|9]we show the SFR/^/SFR^q distribution. In the 
figure it can be seen that Ha estimates are systematically 
above the IR estimates. However, for 91% of the objects, 
the star formation rates agree within a factor of 3. If 
we consider only our spectroscopic confirmed sample we 
find very similar results (see table [1]): (SFR# Q )=17tio 
Moyr" 1 and (SFR ZJl /SFR ffa )=0.96. 

5.3. Exploring the scatter 

It is interesting to explore the reasons that originate 
the discrepancies between different tracers. The Ha line 
is only produced when the star forming region includes 
stars with masses above 10 Mq. Thus, only star forming 
regions aging less than 20 Myr are detectable through 
this line, since older regions would not have stars mas- 
sive enough to photoionize the surrounding gas. There 
are other factors that could affect, such as metallicity, 
fraction of ionizing photons that escape, etc. This set 
of conditions do not hold for the other tracers, which 
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Fig. 8. — Comparison of SFRs inferred from Ha luminosity and 
IR luminosity. Top: No extinction correction applied. Bot- 
tom: Extinction correction applied to Ha. Filled circles are ob- 
jects confirmed by optical spectroscopy, whereas open circles are 
objects without spectroscopic confirmation. We show also the 
SFR/ij/SFRiJQ ratio versus SFRj/a- 

have their own. Ultraviolet, for example, is more sen- 
sible to less massive stars, being more affected by the 
star formation history of the galaxy. Infrared is also 
affected by the star formation history as evolved stars 
can make a significant fraction of the infrared emission 
<jda Cunha et al.ll2008t l. iKennicutt et al.l (|2009f l find that 
50% of the infrared emission in the local galaxies within 
the SIN GS sample come fr om evolved stars. At higher 
redshifts, Sali m et al.l ([2P09) find for a sample with z<1.4 
and with star formation (NUV-R<3.5), that the infrared 
flux fraction originated by intermediate and old stars can 
be as high as 60%, with a typical value ~40%. 

The calibration of the different tracers implies the as- 
sumption of a star formation history. This could lead 
to discrepancies when comparing different tracers. As a 
measure of star formation history we use the Ha equiv- 
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Fig. 9. — Histogram of the ratio SFRrfj/SFRjja in logarithmic 
scale. Most of the objects concentrate around the unity ratio, 
though in general the iR estimates are higher than those obtained 
with Ha. The dashed lines encloses the objects whose SFRs agree 
within a factor of three. 

alent width. It is interesting to check if there is any 
systematic difference depending on age or star forma- 
tion history. The equivalent width is defined as the 
quotient between the Ha flux, which measures the rel- 
evance of the recently formed population, and the con- 
tinuum flux under the line (by wavelength unit), which 
measures the contribution of the stars formed before. 
It is also related to the age of the star forming region 
(|Perez-Gonzalez et al.l [2003) . In figure [TU] we show the 
SFRpfyv/SFR^Q, ratio versus the Ha equivalent width. 
There exists an anti-correlation among both magnitudes, 
with the best fit given by: 

M ^SJf^ ) = (1.45±0.64M0.72±0.29)xlog(W Ha ) 

O-t tin a 

When the equivalent width is low, the weight of the 
young stars is less significant compared to that of the 
old population. In this case, the UV provides higher star 
formation rates than Ha, as it is more sensitive to older 
stars. As we move towards higher equivalent widths, the 
recently formed stars become more and more important 
over the old population. For the higher EW values the 
UV subestimates the SFRs compared to Ha. 

If we repeat this methodology regarding Ha and the 
IR, we obtain similar results. An anti-correlation exists 
between the ratio SFR/^/SFR^q, and EW(Ha) (see fig- 
ure [TTJ, which is described by the best linear fit as: 

log( S JJl IR ) = (1.09±0.76)-(0.52±0.35)xlog(£W/ ffa ) 
bb tin a 

(12) 
Again, as we move towards higher EW(Ha), i.e. to 
higher contributions of young stars, the SFR/#/SFR# Q 
ratio decreases. We have also plotted the local UCM 
sample of star forming galaxies. These galaxies are lo- 
cated in the same region as the z~0.84 sample, so the ef- 
fect of the star formation history on the SFR/ij/SFR/f Q 
ratio is similar at both redshifts. In both cases, con- 
straining to the spectroscopically confirmed sample, 
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Fig. 10.— SFR FC /y/SFRH Q ratio versus EW(Hq). The line 
is the best linear fit to the data in logarithmic scale. The 
high EW(Ho) objects tend to have lower SFR^-jjy/SFR^jQ, ra- 
tios whereas the lower EW(H«) objects have tend to have higher 
ratios. 

yields compatible results within errors. We note that 
the significance of both relations is not very high as er- 
rors are large, due to measurement errors, different star 
formation histories, etc. 

Several authors (iPerez-Gonzalezl 120031 : iFlores et al.l 
120041: lHammer et al.l 120051) found a correlation between 
this ratio and the IR luminosity in the Local Universe. 
The authors argued that the more luminous is a galaxy 
in the infrared the more affected by extinction, to the 
point that the optical tracers could lose an important 
fraction of star formation, with some regions totally 
obscured by extinction. This would explain the sub- 
estimation of the SFR when measured by optical esti- 
mators, even when applying extinc tion corrections. At 
higher redshift ICardiel et al.l (|2003| ) found a similar be- 
havior for their sample at z=0.8. In figure [12] we repre- 
sent the SFR//j/SFR# a ratio versus the infrared lumi- 
nosity L/fl(8-1000/im). We find the same behavior re- 
ported by these authors: the Ha estimator starts to sub- 
estimate the SFR (with respect to IR) when we move to 
higher Lrfl(8-1000/im). The dependency it is not a selec- 
tion effect as the limits of our sample would allow us to 
detect galaxies with higher and lower ratios (see figure). 
No dependency is found between the extinction of the 
objects (coded with different colors in the figure) and 
the degree of sub-estimation. One would expect some 
kind of dependency, as the obscured regions in the op- 
tical are visible in the IR, although it does not discard 
this scenario. 

The best linear fit, excluding the 5% extreme values, 
is given by: 



log (Hf^) = -(2.61±2.17)+(0.23±0.19)xlog(^ 

(13) 

which can be seen in the figure as a thick line. 

This result agrees with previous results, with the IR 
providing higher SFRs as we move towards higher IR lu- 
minosities. We note that uncertainties are large and in 
this case, where we are representing x/y vs. x, a small 
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Fig. 11.— SFRja/SFRhc ratio versus EW(Ha). The line is the 
best linear fit to the data in logarithmic scale. The high EW(Hc?) 
objects tend to have lower SFR/^/SFRj^a ratios whereas the lower 
EW(Ha) objects have tend to have higher ratios. The crosses rep- 
resent the UCM sample of local star forming galaxies, which show 
the same trend. 

correlation would arise from random scatter. However, 
the slope of such correlation is always below 0.1 within 
the range of our values, as we have checked generating 
random samples. We have plotted the UCM local sample 
as a reference, with its best linear fit as a dotted line in 
figure [T2"1 Both samples present very similar slopes when 
fitted, but there is an offset between them. For the UCM 
sample the IR starts to provide higher SFRs for lower 
IR luminosities than for the z~0.84 sample. If we con- 
sider that the change in the SFR/fl/SFR# Q ratio is due 
to the increment of Ha luminosity totally obscured by 
dust, at higher IR luminosity more regions would be to- 
tally obscured by dust. Then, the difference between the 
local relation and that at z^0.84 presented in Figure IT21 
could be explained by a change in the number and size 
of star forming regions. These should be less numerous 
at z=0.84, although larger, given that SFRs are higher 
in the sample at z=0.84 than in the local sample. Thus, 
Ha luminosity would be higher for each region and the 
dust would not totally attenuate that region. Only in 
very luminous galaxies in the IR with large amounts of 
dust it would be possible to totally attenuate some star 
forming regions in the visible. 

Thus, we have two possible effects that could explain 
the scatter when comparing SFRs coming from Ha and 
IR: a) contribution of the evolved population and b) pres- 
ence of star forming regions totally attenuated by dust. 

However, figurc[T2]can be explained as well taking into 
account the contribution of the evolved population to 
the IR. We have shown that the SFR /fi /SFR ffa ratio in- 
creases with the IR luminosity and that there is an offset 
in that relation between the local Universe and z~0.84. 
But, the dependency of that ratio with the EW(Ha), i.e. 
with the weighted age, is independent of redshift. Then, 
we can consider that the same SFR/^/SFRff Q ratio im- 
plies the same weighted age at z=0 and at z^0.84. On 
the other hand, galaxies with the same weighted age have 
higher infrared luminosity at z=0.84 than at z=0, as de- 
picted in figure [T2l As their IR luminosity is higher, the 
star formation rate is higher, and, as the weighted age is 
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Fig. 12. — SFR/b/SFRj^q ratio as a function of IR luminosity. 
The extinction is color coded as can be seen in the legend. The 
filled circles are confirmed by spectroscopy whereas the open circles 
are not. The thick line represents the best linear fit to the data 
(see text). The red line is our selection limit. We could detect only 
objects below this line. The UCM local sample of star forming 
galaxies is represented by the crosses. The dotted gray line is the 
best linear fit to this sample. 

the same, the underlying population has to be more lumi- 
nous, i.e. more massive. This is in ag reement with what 
can be expected from the downsizing (jCowie et al.lll996f) 
scenario, where star formation moves from more massive 
galaxies at higher redshifts to less massive galaxies at 
lower redshifts (see section [S]). 

Although the effect of age could explain the observed 
difference, probably both age and extreme attenuation 
of some regions are contributing. 

6. STELLAR MASSES 

The stellar mass is one of the most important prop- 
erties of a galaxy, as it provides a robust measurement 
of the scale of the galaxy and it is also an indicator of 
the past star formation. The estimate of stellar mass 
is obtained from the best fitting template to the SED 
of each galaxy. The template provides mass-to-light ra- 
tios for each observed band and a stellar mass is com- 
puted for each one. The final value is the average of 
the values obtained for each observed band, being the 
associated error the standard deviation of the distribu- 
tion of stellar masses. The results are more reliable than 
those obtained through a single mass-to-light ratio, as it 
is less sensitive to the star formation history or errors in 
photometr y or templates. For more de tails on the pro- 
cedure see IPerez-Gonzalez et al.l (2008) and IBarro et al.l 

unns). 

Stellar masses were obtained with the PEGASE 2.0 
(|Fioc fc Rocca-Volmerangelll997fl stellar population syn- 
thesis models, a Salpeter (119551) initi al mass function 
(IMF) and the iCalzetti et al.l (I2000D extinction law 
Differ ent stellar popu lation models (jBruzual fc Chariot] 
12003 iMaraston l I2005L Chariot & Bruzual 2009), IMFs 
(jKrou pa 2001; iChabrierl 120031 ) or extinction laws will 
provide different estimates. The models used here pre- 
dict the largest stellar masses, although all models are 
roughly consistent within a factor of 2. For a detailed 
comparison between the different models we refer the 
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Fig. 13. — Histogram of stellar masses for our sample. The com- 
pleteness falls below 10 1() M©. 



reader to Bar ro et all (|2011b[ ). 

In figure [13] we present the histogram of masses for 
our sample. The median and standard deviation for the 
distribution are M ir =lAtf,l x lO^ M^. At the same red- 
shift, the typical mass found by IPerez-G onzalez et al.1 
(2008) for an IRAC selected sample is M*=1.6 x lO n M . 
The typical mass of an star-forming galaxy (to the limit 
of our sample) is ten times lower th an the typical mass 
of the global population of galaxies. iSobral et al.l (J2010D 
(hereafter S10) find a typical mass M*=2.25xl0 10 M Q 
(after scaling from a Chabrier to a Salpeter IMF) for 
their HiZELS sample at z=0.84, in very good agreement 
with ours, given that both limiting fluxes are very simi- 
lar. 

The loss of the low mass population of star forming 
galaxies is clear in the histogram, with the number of 
galaxies starting to decrease below ~ 10 10 M Q . This 
effect is produced by the limiting line flux reached in 
our selection process, given the correlation between stel- 
lar mass and star formati on found in the Local Universe 
(|Brinchmann et al.l 12004). hereafter B04) and at higher 
redsh ifts up to z~6lNoeske et all I2007bt lElbaz et al.l 
2007J; iDaddi et al.l 120071: IStark et alj|2009ft . 

The SFR-M* correlation for our sample is shown in 
figure Q31 where the Ha SFR versus the stellar mass is 
represented. The completeness level (red dashed line) 
is estimated performing simulations of the whole selec- 
tion process and taking into account the extinction. This 
process involves several steps which are explained in de- 
tail in appendix \M T he SFR-M* relations obtained 
from Dut ton et al.l (2010). from the s amples presented in 
INoeske et al.l (|2007bD at z~0.8 and in lElbaz et all (|2007f ) 
at z^l.0, in good agreement with our sample, are over- 
plotted. In a ddition, the UCM local samp le of star form- 
ing galaxies (jPerez-Gonzalez et al.|[2OO30 and the SDSS 
DR4 galaxies classified as star forming in B04 are also 
shown. The slope for these samples are similar to ours, 
although the SDSS is steeper. There is a shift in SFR be- 
tween the local ones and that at z~0.84. For a given mass 
the sample at z^0.84 presents higher (~ x5.5) star for- 
mation rates than the local sample. Another way to see it 
is that for a given star formation rate, the galaxies in the 



past were less massive than local galaxies. This difference 
between the local Universe and z~0.84 clearly shows that 
star formation changes as the Universe evolves. Con- 
trary to our result, S10 do not find any relation between 
SFR and stellar mass. This is intriguing as the selection 
technique and line flux reached are very similar in both 
surveys. 

In this type of comparisons between our sample and 
the local samples one could think that we are just watch- 
ing the tail of the distribution for the z~0.84 sample, 
which lead us to the wrong conc lusion that an evolu- 
tionary effect is present. In IV08I we computed the Ha 
luminosity function and we also estimated the complete- 
ness limit. The results showed that we start to lose a 
substantial fraction of objects one order of magnitude 
below L*(Ha) (obtained from the LF fit to a Schcchtcr 
function). We were 50% complete for log L(Ha) > 42.0, 
which in terms of L*(Ha) is log(L(H a) /L* (Ha)) ~ -1.0. 
Therefore, we conclude that we are not observing the 
rare galaxies in the tail of the distribution. On the other 
hand, the UCM sample is considered complete down to 
log L(Ha) ~ 40.7, which in terms of L*(Ha) (for this 
redshift) is log (L(Ha)/L*(Ha)) ~ -1.2. Thus we are 
reaching very similar Ha luminosities for both samples 
in terms of L*(Ha). Another argument supporting this 
conclusion is that we detect essentially the same popula- 
tion as_samples selected in UV or IR at the same redshift 
(see Ell) 

It is also interesting to compare with a sample 
purely selected in stellar m ass, as that presented in 
IPerez-Gonzalez et al.l ((2008) . The median values of this 
sample at z=0. 8-1.0 fall very close to our best fit, be- 
ing compatible within errors. However, the difference 
between our best fit and the median values of that sam- 
ple increases as we move to higher masses, presenting 
the mass selected sample lower SFRs. Although this 
trend is very weak and, given the errors, could even not 
be present, is con sistent with the downsizing scenario 
(|Cowie et al.l [1996) : the fraction of galaxies with very 
low star formation seem to increase as we move to higher 
masses, decreasing the median SFR in the mass selected 
sample. At low masses the effect is less pronounced as 
galaxies arc undergoing strong star formation episodes 
and are still very active in general. 

The SFR-M* correlation allows us to estimate the mass 
range in which we can consider our sample unbiased. 
The cut between the SFR 50% completeness level and 
the lower envelope in the SFR-M* distribution gives us 
the stellar mass range within we can consider our sam- 
ple free of biases. The lower envelope is the best linear 
fit shifted downwards to enclose 95% of the data con- 
firmed by spectroscopy (90% within lower and upper en- 
velope). In practice it is an iterative process: first we 
fit the data above and below initial mass limits, then we 
find the new mass limits in the intersection between the 
lower envelope and the completeness curve, and the pro- 
cedure is repeated until the mass limits used in the fit 
and the ones obtained converge. This gives us a lower 
limit of ~1O 1O M0 and no upper limit, indicating that we 
are limited by explored volume on the upper side, as we 
do not detect any galaxy above ~3xlO n M0. As men- 
tioned bef ore, the derived corr elati on is in good agree - 
ment with INoeske et al.l (|2007b[ ) and lElbaz etall (|2007l) . 
both at a similar redshift. Moreover, we find that the 
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Fig. 14. — H« star formation rate versus mass. The circles represent our sample: filled for the spectroscopic confirmations and open for 
the remainder. The red dashed line represents the 50% completeness level. The best fit within mass c omp leteness li mits is sh own as a solid 
line, wi th the dashed lines enclosing 90% of the data. The filled region are the best fits obtained by Elbaz et al. (2007) and Noeske et al. 
(2007b) (compiled in Dutton et al. 2010). The mass selected sample values at z=0. 8-1.0 obtained by (Perez-Gonzalez et al. 2008, PG08) 
are represented by the blue crosses. The UCM (crosses) and SDSS (blue colormap) local samples are also shown. 



scatter is ^0.3 dex, which is the typical value found in 
other studies and seems t o be almost constant with red- 
shift (jDutton et al.|[2010l) . However, the lower envelope 
and the 50% completeness limit are close, and thus it is 
still possible that the observed correlation is produced by 
the selection effect. To rule out this possibility we simu- 
late fake samples of galaxies following different SFR-M* 
relations. First we simulate a population of galaxies fol- 
lowing the linear relation derived from our sample. We 
assign to each galaxy a random stellar mass. Given this 
stellar mass we compute the SFR with the linear rela- 
tion derived from the real sample, adding gaussian noise 
to simulate the scatter (<r=0.3). Once we have this fake 
population we check if the galaxies would be detected 
considering the completeness curves (measured for differ- 
ent completeness levels) and the volume sampled. After 
repeating the fitting process considering the complete- 
ness curve as well as the lower envelope, we obtain that 
the results are in good agreement with the input, with 
relative errors within 10% for the slope and 20% for the 
constant term. The same process is repeated 20 times 
to avoid biases due to rare distributions. When a flatter 



relation (20% less steep) is used we find similar results, 
although in some cases the lower envelope is so low that 
it does not cut the completeness curve and no measure 
can be obtained. Thus, given the SFR-M* relation and 
the limits of our sample we can be confident that we will 
not introduce any substantial bias when inferring mass 
related properties using the sample within those mass 
limits. 

Figure [T4l could be interpreted as going against down- 
sizing, as galaxies at z~0.84 have higher SFRs than local 
ones, independently of their stellar mass. The key con- 
cept here is not the absolute SFR, but the sSFR, which 
is the SFR per unit of stellar mass, and thus is a good 
indicator of the impact that the star formation has in the 
galaxy. In figure [15] we represent the sSFR for our sam- 
ple and for the SDSS and UCM local samples. Now it is 
clear the change in star formation as we move towards 
the local Universe, as well as the shift in star formation 
from more massive to less massive galaxies, considering 
the evolutionary impact of the SFR processes. 

There exists an anti-correlation between the sSFR and 
the stellar mass, evidencing that star formation processes 
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have a higher impact on less massive galaxies. This trend 
is also present in the star forming galaxies of the local 
Universe, with a similar slope in the case of the UCM 
sample, although shifted in the sSFR axis to lower val- 
ues, indicating that the star formation is less important. 
Most of our sample fall below the line (green dotted line) 
where the L*(Ha) galaxies would lie at that redshift, 
which shows that we are not missi ng the g eneral p opula- 
tion. Previous determinations (e.g|Pere~G onzalc z et al.l 
12001 INoeske et al.ll2007at IRodighiero et al.ll2010l S10) 
already found this relation but discr epancies arise re- 
gardin g the slope of the correlation. IRodighiero et al.l 
(|2010f ). through IR SFR estimat i ons, fo und a flatter re- 
lation compared to INoeske et al.l (|2007al ). who used UV- 
optical SFR estimators. On the contrary, S10 found a 
much steeper relationship, with a slope ~ -1, given that 
these authors do not find any correlation between SFR 
and mass. 

The slope for our sample (computed only for objects 
within our mass limits ) is /3=-0.4±0.1, whi ch is in good 
agreement with that of lNoeske et al.l (l2007bD . but s teeper 
than the value obtained bvlRodighic ro et al.l (|2010( ) (/3=- 
0.28) in the redshift range 0.5 < z < 1.0. This dis- 
crepancy could come from the different selection criteria. 
While in t his paper galaxies are selected purely by star 
formation, Rodighi ero et all analyzed a sample selected 
in mass through the IRAC 4.5/Ltm band, with a color cut 
to avoid the inclusion of red sequence galaxies, although 
they try to recover dusty starbursts that might fall in 
the red sequence. In this work we have not excluded 
any galaxy by its color since, as we have shown, most of 
the red-sequence galaxies were in fact dus ty star form- 
ing ga laxies. In line with this hypothesis, iKarim et all 
(|2011[ ) found a similar slope to ours (/3=-0.38), using a 
sample of star forming galaxies selected from a mass se- 
lected sample in the IRAC 3.6/im band. However, in 
this case the authors classify the objects as star-forming 
if they belong t o the blue cloud, once the attenuation has 
been removed. Gilbank et al.l (|2011l) have found very re- 
cently a very similar value for the slope (/3—-0A2) using 
a spectroscopic sample with [Oil] SFRs around z~l. 

The birth rate parameter b is linked to the sSFR: 



b = 



SFR 

(SFR) 



= SFR- 



tf 



2xM, 



(14) 



where (SFR) is the average star formation in the whole 
history of the galaxy, or, in other words, the average 
SFR that would have produced the current stellar mass. 
Thus, we can obtain this average star formation divid- 
ing the stellar mass M* by the elapsed time since the 
galaxy formed tf. The factor two takes into account the 
stellar mass returned to the interstellar medium. The 
parameter depends on the choice of the beginning of star 
formation, but it is still very interesting when comparing 
populations at different redshifts. In this work we set 
this initial time at the beginning of the Universe, t/=6.5 
Gyr for z=0.84 and t/=13.4 Gyr for the local Universe. 
A galaxy with a value of b higher than one tell us that 
the current star formation is more intense than the aver- 
age star formation in the past. A value of two indicates 
a specially intense star formation episode. 

Most of our galaxy sample presents values of b higher 
than one, with 85 % (119/140) with b>l and 66 % 
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Fig. 15. — Specific H« star formation rate versus mass. The cir- 
cles represent our sample: filled for the spectroscopic confirmations 
and open for the remainder. Our 50% completeness level is repre- 
sented by the red dashed line. The diagonal dotted line is the place 
that would occupy the objects with L(Ha)=L*(Ha). The crosses 
represent the UCM sample of local star forming galaxies. The color 
map shows the values for the SDSS sample. The horizontal dashed 
lines indicate the values of sSFR at which 6=1 and 6=2 at z=0.84. 
The dotted lines are the same but for the local Universe and 6=1. 
The line for 6=2 coincide with the 6=1 line for z=0.84. 

(92/140) with b>2. In the local Universe the SDSS and 
UCM samples show a very different scenario. Half the 
SDSS sample (48%) have b values over one and only 16% 
above two. The UCM sample presents similar results, 
with 58% of the sample with b>l and 25% with b>2. 
If we confine the analysis to the most massive galaxies 
(with stellar masses above 10 11 M Q ) within our sample, 
36% (5/14) of these have b values above one, and 29% 
(4/14) have values over two. Regarding the local sam- 
ples, the proportions are very different, with 20% in the 
case of the SDSS and 14% in the case of the UCM with 
b over one, and only 3% with b over 2 within the SDSS 
sample and none in the UCM sample. This is a direct 
evidence of downsizing, as the fraction of most massive 
galaxies with intense (&>1) or very intense (b>2) star 
formation was higher in the past, and it has reduced 
dramatically from that epoch to the present. 

6.1. Quenching Mass 

The star formation mass relation holds up until a cer- 
tain mass, above which it no longer holds and the star 
formation drops sharply (see for example B04) . Galaxies 
above this mass are considered quiescent, as star for- 
mation processes are no longer the main drivers of its 
evolution and they move to the red sequence. Therefore, 
this mass defines an upper limit to the stellar mass o f 
the galaxies actively forming stars. iBundv et al.l (2006) 
(hereafter B06) found that this quenching mass evolves 
with redshift, increasing as we move to higher redshift, 
as expected in the downsizing scenario. In this work, 
we use our Ha selected samples to study this quench- 
ing mass evolution, which our star forming samples also 
reflect. The decrease of sSFR with mass implies that 
galaxies will eventually reach a mass over which the star 
formation processes will be very low and we can consider 
them quiescent. The observed shift between both trends 
leads to a different quenching masses, being lower that 
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Fig. 16. — Doubling time t^ versus stellar mass. The circles rep- 
resent our sample: filled for the spectroscopic confirmations and 
open for the remainder. The 50% completeness level is represented 
by the red dashed line. The crosses represent the UCM local sam- 
ple of star forming galaxies. The dashed horizontal line indicates 
the doubling time at which we consider a galaxy at z=0.84 as qui- 
escent. The dotted horizontal line represents the same but for the 
local Universe. Best linear fits computed for both samples are also 
shown. 

at the local Universe. 

With our data it is possible to estimate an upper limit 
for the quenching mass. For the sake of clarity we are 
going to use the doubling time td, which is analogous to 
the sSFR and is defined as: 



t d 



M» 



1 



SFR (1 - R) SSFR (1 - R) 



(15) 



where R is the fraction of mass returned to the inter stel- 
lar m edium which is generally assumed to be ~0.5 (j Belli 
l200l . 

The doubling time tells us how long will take for that 
galaxy to duplicate its stellar mass if its current star for- 
mation stays constant. Galaxies with a large doubling 
time will evolve slowly whereas galaxies with a small 
one will evolve quickly. Doubling times versus mass for 
our sample and the local UCM sample are shown in fig- 
ure llGI In order to estimate the quenching mass we define 
a galaxy as quiescent if its doubling time is higher than 
what we define as quenching time: tq = 3 x i#, where 
tii is the Hubble time. To obtain the typical mass which 
corresponds to the quenching time we performed several 
steps. First, we simulated 1000 realizations of our sam- 
ple, varying randomly the values of SFR and mass within 
twice the errors, i.e., each object will have values ran- 
domly distributed in the intervals [M-2AM, M+2AM] 
and [SFR Htt -2ASFR HQ , SFR Hq +2ASFR Hq ]. Second, 
we do a linear fit of td versus mass only with the objects 
whose simulated mass fall above our mass limit. For each 
of these fits we compute the quenching mass as the mass 
at which the doubling time td is equal to the quenching 
time tQ. The final quenching mass Mq is the median 
of the whole distribution of quenching masses, with the 
error determined by the standard deviation of the dis- 
tribution. The same process has been followed for the 
UCM sample. 
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Fig. 17. — Evolution of the quenching mass limit Mq. Red cir- 
cles represent the results obtained in this work for the Ha selected 
samples at z~0.84 and the local Universe. The orange circles are 
the estimated quenching masses when no extinction correction is 
considered. The rest of the points correspond to the B06 work ac- 
cording to the different criteria employed: black squares for mor- 
phology, green crosses for [OIIJA3727A SFRs and blue triangles for 
the (U-B) color. 

(log (M Q /M )=12.O±O.2) for the z-0.84 sample and 
M Q =7.9tl;5XlO lo M (log(M Q /M )=lO.9±O.l)forthe 
local sample. If we consider only the spectroscopically 
confirmed sample we obtain log (Mq/M )=12.2±O.2, 
slightly higher although compatible within errors. In the 
case of z^0.84 the quenching mass is outside the range 
of masses detected, given the limit on the detection of 
massive galaxies imposed by sampled volume and the 
equivalent width limit of the survey, which prevents us 
from selecting objects with lower sSFRs. These masses 
are upper limits, given that at high stellar masses 
the correlation between doubling time and mass will 
break as a consequence of quenching. Galaxies with 
higher td than predicted by the correlation will appear 
as the quenching takes over, possibly lowering the 
average quenching mass, specially in the case of z~0.84, 
where no galaxies around Mq have been detected. In 
order to detect these galaxies it would be necessary to 
survey larger volumes. In addition, the simulations (see 
section [5]) show that we may overestimate the quenching 
mass at z~0.84 by ~0.1 dex, due to the completeness 
limits. However, these simulations also shows that we 
would be able to detect quenching masses ~0.5 dex lower 
(SFR-M* slope 20% lower), with a similar dispersion. 

Our quenching mass estimation for the local Universe is 
in very good agreement with the stellar mass (~7xl0 10 
M ; scaled from a Kroupa t o a Salpeter IMF) above 
which Kauffm ann et al.l ([2003D found a rapid increase in 
the fraction of galaxies with old population in the SDSS 
local sample. This change, detected by a transition from 
lower values of D„(4000) to higher values, is also seen as 
a change in the slope of the /i*-M, (surface stellar mass 
density versus stellar mass) correlation. Our result at 
z~0.84 is higher than those estimated by B06 at a simi- 
lar redshift. In their work, they used three different ap- 
proaches: morphology, U-B color and SFRs derived from 
[Oil] equivalent width. Through the morphology crite- 
rion they obtained Mq ~8xlO 11 M whereas both color 
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TABLE 1 

Comparison between total and confirmed samples. 



Property 
(1) " 


Total Sample 
(2) 


Confirmed Sample 
(3) 


7 Calzetti et al.(2000) 
{SFR Ha ) 

{SFIWv/SFRhc} 

<SFR /fl /SFR Ha } 

log (Mq/M q ) 


0.55 

lltf MQyr" 1 

0.89 

0.95 

12.0±0.2 


0.56 

14if Moyr" 1 

0.87 

0.96 

12.2±0.2 



Note. — (1) Measured property (2) Value obtained using the 
whole sample (3) Value obtained using the sample confirmed 
with optical spectroscopy 

and SFR criteria provide lower masses ^lO 11 M (see fig- 
ure [T7]). We have scaled B06 masses from the Chabrier 
IMF to the Salpeter IMF used in this work adding 0.25 
dex. Our result is in good agreement with their morphol- 
ogy based estimation but it is higher than those based in 
color or SFR. B06 attributed this difference to a longer 
timescale in the processes that transform late types into 
early types. However our value is solely based in star 
formation and no morphology considerations have been 
done. One of the caveats of their SFR and color measure- 
ments is that extinction was not corrected. Therefore, 
dusty starbursts would appear redder and with lower 
SFRs, as they would be classified as red or non-star form- 
ing galaxies, which translates in lower quenching masses. 
If we estimate again Mq for our Ha selected sam- 
ples, but this time without applying the extinction cor- 
rection, we obtain lower values: MQ—7.6~til : xW 10 M ( 7 ) 
log M Q /M Q =10.88±0.09) for the z-0.84 sample and 
M Q =1.3+^2Xl0 10 M Q (log Mq/M q =10.13±0.05) for 
the local sample. The effect of the extinction is very 
high, and is enough to account for the difference between 
BOG morphology and color/SFR e stimations. Our result 
is also consistent with the work bv lPerez- Gonzalcz et al.l 
(2008), where they found that at z~0.8 all the stellar 
mass has been already assembled for objects more mas- 
sive than M*=10 12 M Q and almost fully assembled for 
objects with stellar mass in the range 1O 117 -1O 12 M0. 

The quenching masses estimated through this method 
rely on the definition of tQ, however, independently of 
this parameter, we find a strong evolution between the 
local Universe and z~0.84. In the local Universe, galaxies 
with mass higher than 10 11 M Q are quiescent and their 
evolution is limited to interactions with other galaxies 
via dry mergers (or other processes not involving massive 
star formation), whereas at z^0.84, galaxies with mass 
in the range ^lO 11 — ^1O 12 M are still under strong star 
formation processes. This is, again, in good agreement 
with the downsizing scenario. Despite we do not have 
enough data to constrain the Mq evolution with red- 
shift we find that our results are compatible with the 
parametrization given by B06, i.e Mq oc(l+z) 45 . The 
added value is that we have extended the redshift base- 
line to the local Universe, using samples selected uni- 
formly. 

7. SUMMARY AND CONCLUSIONS 

In this work we have analyzed the properties of an 
Ha selected sample of star-forming galaxies at z~0.84, 
focusing on the star formation and stellar mass. 



We have discarded the AGN contaminants through 
two criteria: X-ray luminosities and IRAC colors. We 
find seven counterparts in X-ray, though three of them 
present very low fluxes, compatible with being originated 
by star formation processes. Thus, we only discard the 
four objects with fluxes high enough to have an AGN ori- 
gin. Using IRAC colors we find another ten objects (one 
of them already detected in X-rays) that fulfill the cri- 
terion to be considered AGN. A total of thirteen objects 
are finally discarded. 

The objects of our sample present a median Mb=- 
20.5±0.9 m , brighter by more than one magnitude than 
the UCM local sample of star forming galaxies. Most 
of the galaxies belong to the blue sequence with a small 
fraction of objects in the green valley and the red se- 
quence. Once the extinction corrections are applied all 
these red objects except two move to the blue sequence, 
unveiling their dusty nature. 

A check o n the extinction law reveals that the 
iCalzetti et al.l (|2000D extinction law is appropriate for 
our objects, but with E(B-V) stars =0.55 E(B-V) ffas . 

The Ha SFR, without applying the extinction correc- 
tion, presents values in the range 2-10 M Q . We have also 
estimated SFRs with FUV and IR. In the first case, the 
non-extinction corrected FUV underestimates the SFR 
with respect to Ha. The opposite case is given for the 
IR, which overestimates the SFR with respect to Ha. 
These discrepancies are mainly driven by the extinction. 
Once we apply the extinction correction to both FUV 
and Ha estimations, all SFR tracers agree within a fac- 
tor of three and the highest SFRs reach several hundreds 
solar masses per year. 

The scattering between the different tracers are cor- 
related with the Ha equivalent width. This can be ex- 
plained through the different weighted age of the objects 
(which is related to the EW) and the fact that FUV an 
IR SFRs are sensitive to a longer time range than Ha, 
being more affected by older populations. 

We have estimated stellar masses for our sample, find- 
ing that the median value is M + =1.4+^ x 1O 1O M , in 
good agreement with the result obtained by S10 for a 
sample selecte d with similar criteria as ours . The typical 
mass found by Perez-Gonz alez et al.l (|2008f) for an IRAC 
selected sample at this redshift is ten times higher than 
this value. 

Our sample shows a trend between SFR and stellar 
mass. The slope of thi s trend is in good agr eement with 
the value o btained bvlNoeske et all (|2007a[) and is flat- 
ter than the lElbaz et al.l (|2007l ) result. The trend is very 
similar to that of the local Universe, although shifted to 
higher values of SFR. This indicates that, for the same 
stellar mass M*, star forming galaxies at z~0.84 are un- 
der stronger star formation episodes than their local anal- 
ogous. 

The sSFR shows a negative correlation with stellar 
mass. The star formation in more massive galaxies, al- 
though with higher SFR, has less impact than in less 
massive ones, due to the large stellar mass already 
formed. The same trend is observed in the local Uni- 
verse, though shifted to lower sSFRs. This is in good 
agreement with the downsizing scenario, in which mas- 
sive galaxies are formed earlier than less massive ones. 
The fraction of massive galaxies (M* >10 11 M Q ) under- 
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going strong star formation processes (b>2), ^29% at 
z~0.84 against <3% at the local Universe, also supports 
this scenario. 

Finally, we have quantified the downsizing estimating 
the quenching mass at z~0.84 and at the local Universe 
based on the Ha star formation rate. We find that 
M Q =1.0to;4xl0 12 M Q (log (M Q /M )=12.O±O.2) 

at z-0.84 and M Q =7.9+^xl0 10 M Q (log 
(Mq/M )=1O.9±O.1) in the local Universe. The 
evolution since the local Universe is out of doubt, 
with an increase in the quenching mass of an order of 
magnitude. 
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Fig. 18. — Stellar mass versus J band magnitude. The circles represent our sam ple: filled for the sp ectroscopic confirmations and open 
for the remainder. The small dots are the objects within our redshift range from the Barro ct al. (2011a) mass selected sample. The dashed 
line is the best linear fit for the mass selected sample. 

APPENDIX 

COMPLETENESS LEVEL AS A FUNCTION OF STELLAR MASS 

The SFR completeness level for our sample depends on the stellar mass, due to the selection process, which is 
basically a selection in equivalent width. In V08 we simulated the selection process to determine the completeness 
level versus the line flux, but without considering the dependence on stellar masses, as it was unnecessary. The process 
consisted in the introduction of simulated galaxies in the real images to check wether or not they were recovered by 
the selection method. Fake galaxies were created with different line and continuum fluxes, sizes and inclinations in the 
real images. The outcome was the line flux completeness level. 

The problem we need to address here is at which SFR, corrected for attenuation, the completeness level is 50% as a 
function of stellar mass. The main problem is that we need to relate the stellar mass and extinction to the observables, 
i.e. the broad-band and narrow-band magnitudes. Intuitively, one may think that exists a correlation between the 
stellar mass and the broad band magnitude. Indeed, as can be seen in figure [TSl there is a correlation between the 
stellar mass and the J magnitude for our sample. However, it might be affected by selection biases, so we also check 
this relation with another samp le. F or this doub l e-check we use the mass selected sample from the AEGIS database 
presented in lBarro et al.l (|2011al ) and lBarro et al.l (|2011bl ). If we take all the objects with spectroscopic redshift within 
the limits of our sample we obtain a very similar relation. The best fit to the data gives: M*=23.36 - 0.636 mj. 

With this correlation we can assign a stellar mass for an object given its J band flux. Thus, we can check the fraction 
of objects selected for a certain line-flux, defined by their emission in the broad and narrow band filters, and a stellar 
mass, estimated from the broad band flux. The line flux is transformed to SFR applying the corresponding calibration 
(equation [T]) and applying a mean correction for the nitrogen contribution (I([N II] A6584)// (H a)=0. 26). At the end 
we have a completeness level for each combination of SFR and stellar mass. From this we can obtain the SFR 50% 
completeness level as a function of mass. 

However, this completeness does not take into account the attenuation affecting the SFR measurement. To correct 
this effect we can apply the extinction correction to the SFR completeness level: 



C IP ryl A /f *\corrected c IP Dl T\.r*\uncorrected -i n0.4 A(Ha) 

btH{M ) 50% =btK{M j 50% 1U 



(Al) 



uncorrected j s fa e jj a gpj^ f or wmcn we are 59% complete at stellar mass M 
;h 



where SFR(M 
The question that arises now is which amount of extinction to apply. 



The immediate solution is to apply the 
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mean(median) extinction correction obtained for the sample. However, it is well kn own that t he amount o f extin ction 
depends on the total SFR. This fact has been shown for the iSobral et al.l (|2009fl sample by Gam et al.l (|2010l ). In 
particular these authors find that A(Ha)= 0.73 + 0.44 log SFR(IR). As we saw in section [5. 2.21 we can assume that 
the IR derived SFR is equal to the extinction corrected Ha SFR, and thus we can write: 

SFR{M%° Q % ected = SFR{M*)l^° rrected ioO-i(a+bxio g SFR(M'r 5 ^r^ d ) ( A2 ) 

= SFR{M*)l^° rrected ioa'+b'xi°sSFR(M*)r„vr ctad (A3) 

with a=0.73, b=0.44, a'=0.4xa and b'=0.4xb. We can find the value of SFR obtaining: 

l^rv C J? T>( l\J*\uncorrected i „/ 

log SFR(Mr 5 °oT cted = V-y (A4) 

This equation gives us the completeness level at 50% considering that the extinction depends on the total SFR. The 
computed completeness curve is shown in figure [T4l 



